Pattern Formation in a Two-Dimensional Array of Oscillators with Phase-Shifted Coupling 
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We investigate the dynamics of a two-dimensional array of oscillators with phase-shifted coupling. 
Each oscillator is allowed to interact with its neighbors within a finite radius. The system exhibits 
various patterns including squarelike pinwheels, (anti)spirals with phase-randomized cores, and an- 
tiferro patterns embedded in (anti)spirals. We consider the symmetry properties of the system to 
explain the observed behaviors, and estimate the wavelengths of the patterns by linear analysis. 
Finally, we point out the implications of our work for biological neural networks. 
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Large systems of interacting oscillators have been used 
to explain the cooperative behaviors of numerous physi- 
cal, chemical, and biological systems ^2e^M- When 
the coupling between oscillators is sufficiently weak, we 
can describe the dynamics of the system by phase vari- 
ables defined on the limit cycles |^ yj. In the phase- 
oscillator models, the coupling proportional to the sine 
of the phase difference between oscillators has been ex- 
ploited due to the mathematical tractability y, jj]. In 
spite of several successes in explaining the synchroniza- 
tion phenomena, the simple sinusoidal coupling fails to 
account for the collective frequency higher than natu- 
ral frequency 0, dephasing phenomena Q, and effects 
of time-delayed interactions 7]. To resolve these prob- 
lems, phase shifts in the sinusoidal coupling have been 
considered. Most importantly, a nonzero phase shift is 
naturally contributed by the broken odd symmetry of a 
coupling function 8J. The previous works, however, stud- 
ied only the cases of phase shifts in the limited range. 
Moreover, they considered only limited interactions via 
nearest- neighbor coupHng [sIlZIoJ or all-to-all global cou- 
pling |3| which seem to be too restrictive. In the neuro- 
biological systems, for example, it is believed that actual 
coupling takes the forms between these two extremes jlOf . 

In this Rapid Communication, we investigate the ef- 
fect of phase-shifted coupling on the dynamics of a two- 
dimensional array of coupled oscillators. Here we intro- 
duce a finite interaction radius as realistic coupling and 
study over the whole range of phase shifts. We find that 
various spatial patterns come to emerge, and unravel that 
the symmetry properties of the system play an important 
role on the formation of patterns. 

We start with the equations of two coupled phase os- 
cillators 13 

-^ =uj + K sin(6l2 -01 -a) 
at 

^=uj + Ksmiei-e2-a), {lu, K > 0), (1) 

where 6i and 62 represent the phases of oscillators, 
respectively, uj the natural frequency, K the coupling 
strength, and a the phase shift. Phase shift \a\ < 7r/2 



leads to inphasing of the two oscillators, whereas \a\ > 
tt/2 leads to their antiphasing. We find that this sepa- 
ration by |q;| = 7r/2 still holds for an array of oscillators, 
but in a rather sophisticated manner as shown below. 

To investigate the phase-shift effect on the spatially ex- 
tended systems, we study the following model equations: 
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where 6ij denotes the phase of the oscillator at posi- 
tion {i,j) on a two-dimensional lattice, and X^mn = 
Z]mn,o<r„„,j<fl' whcrc rmn,ij IS the distance between 
two oscillators located at (j,j) and {m,n). Each oscil- 
lator interacts with N{R) neighboring oscillators within 
a finite distance R. In the previous works, the nearest- 
neighbor interactions {R = 1) were considered for the 
limited range of a [a, la|- Motivated by nonlocal con- 
nections of neural systems [li|, we explore the general 
cases with i? > 1 as well as i? = 1, over the whole 
range of a. The nonlocal interactions may arise effec- 
tively in the reaction-diffusion systems where chemical 
components constituting the local oscillators are free of 
diffusion while the system involves an extra diffusive com- 
ponent [ijj. 

To reduce complexity of Eq. Q , we perform a trans- 



formation 9ij 
equations: 

de,. 



ujt+Qij, t 



T / K, and get the following 
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y^ sin(e„i„ - Qij - a). 



(3) 



Qij maintains the same spatial patterns as 9ij, but not 
the temporal behaviors (e.g., the phase velocity becomes 
different). Equation (PJ involves only two control param- 
eters, a and R. We numerically integrate Eq. © on a 
rectangular array of 100 x 100 sites with periodic bound- 
ary conditions. We select initial O^ randomly from the 
range [— tt, tt]. 

We investigate emerging patterns as a varies in the 
range [— 7r,7r] for given R (see Figs.^jS)). Based on our 
observations, we separate the cases into those of i? = 1 
and of i? > 1. 




FIG. 1: Phase diagram as a function of a for fixed R. (a) 
The nearest-neighbor coupling {R — 1). al '= 0.297r. S: spi- 
rals, AS: antispirals, T: turbulence, AF-T: antiferro pattern 
embedded in turbulence, AF-AS (AF-S): antiferro pattern 
embedded in antispirals (spirals), (b) The finite-long range 
coupling {R = 6). al '= O.IStt, a2 '= 0.367r, q3 '= 0.577r, 
a4 ■= 0.747r. P: planar oscillation, S-PRC (AS-PRC): spi- 
rals (antispirals) with phase-randomized core, IR: irregular 
pattern induced by expansion of phase-randomized core, NI: 
nearly-incoherent pattern with a weak correlation of a length 
scale ~ R, CPW: competing plane waves occupying their re- 
spective domains of evolvable sizes, PW: plane waves. It is 
also observed that squarelike pinwheels exist transiently with 
plane waves. 



Case R—1. If oscillators are coupled only with their 
nearest neighbors, we get the phase diagram Fig. ^a), 
which shows the symmetric relations between a > and 
a < cases, and between \a\ < tt/2 and \a\ > 7r/2 cases. 

At a = 0, there appear vortices with equiphase lines of 
zero curvature 8]. As a increases, equiphase lines around 
vortices become twisted, and the vortices begin to show 
meandering behavior 13 [Fig.|2fa)]. Further increasing a 
creates many vortex-antivortex pairs and eventually in- 
duces turbulent patterns until a < 7r/2 |5|,ll2| [Fig-I21b)]. 
If we decrease a from 0, similar behaviors are observed, 
but the phase gradients near the vortices become reversed 
such that phases increase radially outward from the vor- 
tices [Fig. |2c)]. In a weak coupling regime K < uj, 9ij 
in Eq. |2Jl grows as time elapses, thus, waves propagate 
inwardly with a < 0. These waves are recently termed 
as "antispirals" [3, llSj . We can understand the antispi- 
ral formation for a < using a symmetry transformation 
a — » —a, Qij -^ — ©ij which leaves Eq. (PJ invariant. In 
other words, the patterns for —a are equivalent to those 



-8,:, for 



Since spirals emerge 



-Qij reverses their phase gradients, 



obtained from O^ 
for a > and O^ 
we get antispirals for a < 0. We also find that the dis- 
persion relation for a < satisfies the known condition 
for antispirals, i.e., the phase velocity is opposite to the 
group velocity Il3l Il4j . 

When R = 1 and |a| > 7r/2, antiferro patterns where 
adjacent oscillators have a phase difference of tt are devel- 
oped, while embedded in the (anti)spirals or turbulence 
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FIG. 2: Patterns obtained from simulations with R—1: (a) 
Spirals (a = 0.27r), (b) turbulence (a — OAn), (c) antispirals 
{a — — 0.27r), (d) antiferro pattern embedded in antispirals 
{a — O.Stt), (e) antispirals by removing antiferro phases from 
(d), (f) the color code of the phase used for all figures. Insets 
of (a), (c), and (d): Magnification of the boxed areas. Arrows 
in (a) and (c) denote the propagation direction of waves with 
9ij in Eq. J^J when K < cu. 



observed a,t a' = a ~tt [Figs.|2Id) and (e)]. An antiferro 
pattern itself, is expected from antiphasing of two cou- 
pled oscillators at \a\ > 7r/2 in Eq. Q. Nonetheless, the 
existence of embedding patterns for a' is rather surpris- 
ing, and we can explain this as follows. By separating Qij 
in Eq. Q into the antiferro phase Qfj^ and the remain- 



ing phase O 
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Consequently, Of^*^ evolves in the same way as Qij for 

a' ~ a — IT in Eq. (Q , and forms the embedding pattern. 

Case R> 1. When R is larger than 1, the nonlocality 





FIG. 4: Patterns realized in a two-dimensional array of the 
Morris-Lecar systems. For the Morris-Lecar equations and 
the parameters, see Ref. |y|. Here we choose / — 0.081, 
k — 0.01. (a) When R — 1, antiferro patterns embedded 
in antispirals are observed. The figure shows only antispi- 
rals after removing the antiferro phases, (b) If R increases 
to 6, coexistence of plane waves and squarelike pinwheels is 
observed. 



FIG. 3: Patterns with R — 6: (a) Spirals with phase- 
randomized cores (a = 0.27r), (b) antispirals with phase- 
randomized cores {a = — 0.27r), (c) irregular pattern induced 
by the expansion of the phase-randomized cores (a = OAn), 
(d) coexistence of plane waves and squarelike pinwheels {a — 
n). 



of interactions imposes several changes on patterns [see 
Figs.HIb) and El . 

First, for \a\ < 7r/2, nonlocal interactions develop 
(anti)spirals with phase-randomized core, contrary to the 
well-defined singularity with R= 1 [Figs. O^a) and (b)]. 
As pointed out in Ref. [ll|, we can understand the in- 
coherence in the core as a result of long-range coupling 
of oscillators which have a broad frequency range due to 
the gradient of effective frequency near the core. 

Second, for |a| > 7r/2, plane waves emerge instead of 
antiferro patterns observed with R = 1. Squarelike pin- 
wheels, where singularities are arranged on a square lat- 
tice, also exist transiently with plane waves [Fig. |3Jd)]. 
These patterns are essentially the same as those observed 
in the models of visual map formation |l5| and time- 
delayed interactions [16| . We observe that these patterns 
have some discreteness for small R > 1 and are clearly 
developed for R > 4. Further increasing R just scales 
the spatial length unit in the continuum limit. We can 
qualitatively explain the appearance of coherent patterns 
such as plane waves. When R > 1, closely located oscil- 
lators within about a distance R share the large por- 
tion of neighbors commonly, thus they experience sim- 
ilar influence. This induces the synchronous movement 
of the closely located oscillators even for dephasing pair- 
wise interactions with |a| > 7r/2. From this viewpoint, in 
the systems such as models of visual map formation [l^ , 
positive short-range coupling in the presence of negative 
long-range coupling may not be crucial for the formation 



of coherent patterns. 

As a validation of our results in a realistic system. 
Fig. ^ displays the results with an array of the Morris- 
Lecar neural oscillators. We take this system because it 
is well known that two diffusively coupled Morris-Lecar 
oscillators exhibit antiphase oscillations as in the case of 
|a| > 7r/2 in Eq. (Q. Actually, the system shows anti- 
ferro patterns embedded in antispirals {R ~ 1) and plane 
waves with squarelike pinwheels {R > 1), as expected 
[Fig. 21 see Figs.j^^e) andj^Jd) for comparison]. 

Now, we want to show that the symmetry argument 
also provides a valuable insight into the quantitative fea- 
tures of the patterns, especially the wavelengths. The 
transformation a —^ a — n, t —f — t leaves Eq. (jSJ in- 
variant, thus the patterns for a' = a — tt are equiva- 
lent to those obtained by t ^ — t for a. In addition. 



-a is equivalent to Qi 



-&ij as already men- 



tioned. Qij -^ ^0jj changes only the direction of a wave 
vector, and does not affect the stability of the mode in 
the isotropic system. Consequently, by the time rever- 
sal, the sign of the growth rate of a linear mode for a 
becomes reversed for a' — ±{a — tt). For |a| > 7r/2, the 
symmetries select only the modes which are excluded for 
al < 7r/2, and vice versa. 

We perform a linear stability analysis to determine the 
stability of the solution of Eq. (Q, ©y (''') = k ■ fij + fir, 
where k denotes a wave vector, n the corresponding fre- 
quency, and fij — [i,j)- We add te^^'^*^^^'^ with suf- 
ficiently small e to the solution and substitute it into 
Eq. ^. The sign of the real part of rj determines the 
growth rate of the perturbing modulation with a wave 
vector ^. If the sign of Re['q) is negative, the perturba- 
tion decays and the solution is linearly stable. By the 
first order expansion and the continuum approximation. 
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FIG. 5: (a) Re{r]) in Eq. © as a function of p. (b) Wave- 
lengths of plane waves (A) as a function of R. Predicted val- 
ues (27r7?/5.1356 — 1.22357?) fit the observed ones well in our 
model for a — tt and in the Moris-Lecar systems for the same 
parameters as in Fig. 2] within relative errors of 0.02 and 
0.03, respectively. The wavelengths observed in our model 
are insensitive to a for the plane-wave regime. 



[H Qi y> IM ll^ ll^ ll^ I and gives an unified perspective 
of pattern behaviors. In general, the presence of phase 
shift in couphng is naturally contributed by the broken 
odd symmetry of coupling functions, and thus our results 
may be applicable in the universal manner. 

We expect that in neurobiological systems, excitory 
(inhibitory) interactions correspond to the case of phase 
shift a > (a < 0), whereas time delay in synapses 
determines whether |a| < tt/2 or not 18]. The specific 
realization of patterns in the systems is left for further 
study. 
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2 cos a 
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where p = Rk and ^ = i?^. We replace sin^[(.^ •r')/2] by 
1/2, the statistical average taken over S's, then obtain 



Reiv) ~ {Re{r]))s 



-Ji(p)cosa, {p=\p\), (6) 



where Ji is the first-order Bessel function. It is verified 
that Eq. © reflects the important properties of Eq. (O 
qualitatively |17l |. 

Figure inia) shows Re{rj) in Eq. ||HJ) as a function of p. 
The sign of Re{ri) for given p is reversed over the a = 
±7r/2 axis, as discussed with the time-reversal argument. 
Re{T]) has the most negative value at p = for \a\ < 7r/2, 
and p '= 5.1356 for \a\ > 7r/2. Therefore, when |a| < 7r/2, 
planar solutions with p = and patterns with p ~ 
are linearly stable, which is consistent with the observed 
wavelengths comparable to the array size. When |a| > 
7r/2, patterns with a wave vector p '=;. 5.1356 are stable 
instead of planar solutions and patterns with p ~ 0. We 
measure the actual wavelengths of the plane waves both 
in our model and in the Morris-Lecar systems, and the 
measured wavelengths turn out to be quite close to our 
predictions [Fig. [S^b)]. 

From the observation that wavelengths with \a\ > 7r/2 
are comparable to the interaction radius R for a two- 
dimensional array, we speculate that dynamics with |a| > 
7r/2 on general interaction networks might reflect the sub- 
structure charateristics of the networks well. 

In conclusion, we have investigated the effect of phase- 
shifted interactions in coupled oscillators in two dimen- 
sions, and found that symmetry properties of the system 
are responsible for the various pattern formation. Our 
simple model succeeds in describing the essential features 
of patterns observed in realistic models and experiments 
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